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Abstract 

There are many different models — both continuous and discrete — 
used to describe gene mutation fixation. In particular, the Moran 
process, the Kimura equation and the replicator dynamics are all well 
known models, that might lead to different conclusions. We present a 
discussion of a unified framework to embrace all these models, in the 
large population regime. 

1 Introduction 

Real world models need to cover a large range of scales. However, models that 
are valid in such a large range are hard to obtain and can be very complex 
to analyze. Alternatively, we might use models that focus on certain scales. 
Thus, on one hand, we might have a microscopic discrete model that is 
derived from first principles while, on the other hand, we might also have 
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continuous models that are easier to analyze but are more phenomenological 
in nature. 

When dealing with many descriptions of the same reality, the connection 
between these various possible descriptions is an important problem. These 
connections, in the particular case of evolutionary game theory for large 
populations, will be discussed here after the work in [7j. In particular, we 
will present a unified theory that covers the Moran process 16j, the mean- 
field theory description by Kimura ^Uj and the replicator dynamics 0. 

In order to develop such a theory, we shall proceed as follows: 

1. we prepare a detailed discrete model — the the MICROSCOPIC model; 

2. we identify suitably scalings and the corresponding negligible (small) 
parameters; 

3. we formally find a new model where these variables are set to zero — the 
THERMODYNAMICAL LIMIT; 

4. we prove that the new model is a good approximation for the previous 
model, within the given scalings; 

5. we study the behaviour of distinguished limiting cases. 

It is important to stress that, usually, step three is obtained from phenomeno- 
logical framework. Thus, even formal connections between discrete and con- 
tinuous models can be very important in understanding their relationship. 
Moreover, this allows one to solve the simpler continuous model and thus 
have the approximate behavior of the solution of the detailed model. This 
approach is classical in the physical sciences where, for instance, continuous 
mechanics can be seen as a formal limit of particle dynamics — although the 
phenomenological derivation has been obtained much earlier pTIj . 

In these derivations, the existence of small parameters is generally natu- 
ral, but the appropriate scalings are not. For example, models of dilute gas 
given by Boltzmann equation converge to the Navier-Stokes or Euler equa- 
tions in fluid dynamics (depending on the precise scaling given) when the 
re-scaled free mean path is set equal to zero [H |2l CH [121 • ^ similar ap- 
proach uses kinetic models for modeling cell movement induced by chemicals 
(chemotaxis) and when the cell free mean path is negligible, their solutions 
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is comparable to the solutions of the Keller-Segel model IH IHl EI- In a 
different framework, relativistic models for particle motion have as the non- 
relativistic limit (i.e., the limit when typical velocities are small compared to 
the velocity of light) the Newtonian physics [131 CS], where quantum equa- 
tions converge again to classical physics when (re-scaled) Planck constant is 
very small [21 E| • 

For the Moran Process, it has been recently noticed that the inverse of 
the population size is the relevant small parameter; cf. H] for instance. 

The outline of this work is as follows: in section |21 we discuss the gener- 
alized Moran process. This includes the standard Moran process as a special 
case, but addresses also the frequency dependent case. In section |21 we re- 
view the scalings and thermodynamical limits found in [7]. The connection 
of some of the thermodynamical limits with the Kimura model is discussed 
in section I3J After that, we briefly outline some of the mathematical issues 
involved in the derivation of the thermodynamical models. Relationship be- 
tween these limits and the Replicator Dynamics is discussed in section IHl 
We then present a series of numerical simulation to illustrate the theory 
discussed and compare results in section [7| Some remarks in more general 
games, where mixed strategies are allowed are given in section |H1 



2 The generalized Moran process 




Figure 1: The Moran process: from a two-types population (a) we chose one 
at random to kill (b) and a second to copy an paste in the place left by the 
first, this time proportional to the fitness. 
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Consider a population of fixed size N, given by two types of individuals, 
I and II. The Moran process is defined by three steps: 

• we choose one of the individuals at random to be eliminated; 

• all the remaining individuals play a game given by the pay-off matrix 



and the individual fitness is identified with the average pay-off. We 
assume that C > A > and B > D > 0. This is the only structure 
of the pay-off matrix that guarantee the existence of non-trivial stable 
equilibrium. This is know in the literature as the "Hawk-and-Dove" 
game. 

• we choose one of the individuals to by copied with probabilities pro- 
portional to the pay-off. 

We repeat these steps until a final state is presumably reached. Intuitively, 
after a long enough time, all individuals will be descendant of a single indi- 
vidual living at time t — 0. More precisely, let P{t, n, N) be the probability 
that at time t we have n type I individuals in a total population of size N, and 
let C-{n,N) {co{n,N), c+{n,N)) be the transition probability associated to 
n ^ n — 1 {n ^ n and n — > 1, respectively). Then, the discrete evolution 
process is given by 



I II 



I 

II 



A B 
C D 



P{t + At, n, N) 



c+{n-l,N)P{t,n-l,N) ( 
+co(n, N)P{t, n, N) + c_(n + 1, N)P{t, n + l,N) . 



(1) 



Let us introduce the vector 



P{t, N) = {P{t, 0, N), P{t, 1, TV), ... , P{t, N, N)y. 



Then the iteration can be written in matrix form as 



P{t + At,N)^MP{t,N), 
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where M is a column-stochastic,tridiagonal matrix. It is also possible to 
show — cf. — that 1 is an eigenvalue of M with multiplicity two, with 
associated eigenvectors given by ei and gat+i. Since the spectrum and its 
multiplicity is unchanged if M is replaced by there are two vectors that 
are kept invariant by M^. One of them is easily seen to be the vector 1 = 
(1, 1, . . . , 1, 1)^. Let F denote the remaining one. Then, [7] showed that F 
yields the stationary fixation probability and also that the quantities 

r/i = (l,P(t,iV)) and r^s = (F, P(t, iV)) 

are invariants of the Moran process. The former is well known and it cor- 
responds to the conservation of probability. The latter, however, seems to 
be new and it can be interpreted as stating that the correlation coefficient 
between a possible state of the Moran process and the stationary fixation 
probability must always be the same of the initial condition. 
We can now prove that 



hm M'^ 
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— Fn^i 




Fn-1 



0\ 
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As a direct consequence, for any normalized initial condition 

Po = (P(0, 0, N), P(0, 1, iV), . . . , P(0, N, N)), 
the final state is 



lim P(KAt, A^) = lim M'^Pq = {1 - A, 0, 0, 



where 



N 



A = J2FnP{0,n,N) 



n=0 



is the fixation probability associated to the initial condition Pq. Note that, 
with probability 1, one of the types will be fixed. This means, that in the 
long range, every mutation will be either fixed or lost. 
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3 Scaling and thermodynamical limits 

The central idea of this section is to find a continuous model that works as 
a good approximation of the Moran process, when the total population is 
large. This means that we want to find a continuous model for the fraction 
of mutants in the limit oo. The core of this process is to define a correct 
scaling for the time-step and for the pay-offs. We will show, however, that 
different scalings will give different thermodynamical limits. But only one 
of these scalings will be able to capture one essential feature of the discrete 
process discussed in the previous section: that genes are always fixed or lost. 
In the continuous setting, this means that, as time increases, the probability 
distribution should move (diffuse) to the boundaries. 

Let us suppose that (formally) there exists a probabihty density 

P(t xN N) 

p(t,x) = lim ^ ' ' ^ = lim NP(t,xN,N) , 
where x = n/N. 

Let us also suppose that this function p : R_,_ x [0, 1] R is sufficiently 
smooth that we can expand the evolution equation as to obtain 

p{t + At,x) -p{t,x) 1 



1 



At NAt 
1 

2 



. 1 



(cj) + + c«) p - (cf - cL°)) d^p] (2) 
+ c^o^ + cL^^) p - (c« - c«) d.p + \ [of + c^^) dip 



ATS 

where ci^ = c*\x), * = +,0, —, i = 0, 1, 2, are defined by 
c+{xN-l,N) = c+(n-l,iV) = cf + lc« + ^cf + 

co{xN,N) = co(n,iV) = cr + ^c« + ^4^^ + o(^ 

(0) , 1 ^(1) , 1 (2) , 1 \ 



c.ixN+l,N) = c4n-l,N) = c'^> + -c'y + ^c'^> + Oi^—j 

So far, we made no assumptions on the behavior of the payoffs A, B, C, and 
D as N ^ oo. Since the appropriate large limit will be attained through a 
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rescaling in time, the fitness should also rescaled accordingly to preserve the 
expected amount of offsprings. In particular, we must have that the fitness 
approaches one as N grows large; see discussion in [7j. Equivalently, payoffs 
must also approach one. The order with they approach unity, however, is a 
free parameter at this point. 
Thus, we impose that 

lim(A5,C,D) = (l,l,l,l), (3) 
lim N%A-1,B -1,C-1, D-l) = {a, b,c,d) , z/ > , (4) 

we find, after a long computation, that 

lim N" (c^^ + + cL'^) = - ^x\a -b-c + d)- 



2x{a ~ c - 2{b - d)) + {d - b) 



lim N' 

N^oo 



- (^cf - cL°^) =a;(l - x){x{a - c) + {1 - x){b - d)) . 



The only non-trivial balances, as can be readily observed, are given by time- 
steps of order At = A^~('^+^), for u e (0, 1]. In this case, we have that 

dtp = -d,{x{l-x){xa + {l-x)P)p) , z/G(0,l) (5) 

and 

dtp = dl {x{l — x)p) — {x{l — x){xa + (1 — x)[3)p) , if z/ = 1 , (6) 

where a = a — c < and f3 = b — d > 0. In the particular case where a = f3 
(i.e., when the fitness is independent of the particular composition of the 
population), the last equation can be shown to be equivalent to a celebrated 
equation in population genetics known as the Kimura equation |TUj . 

The equations above are supplemented by the following boundary condi- 
tions 

1 

p{t,x)dx = 0, for equations (0) and (jH)); 



dt Jo 



d 

4'{x)p{t, x)dx = for equation 



dt 
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In the latter condition, ip satisfies 

tfj" + (13 + {a - (3)x)^' = 0, V'(0)=0 and ^(1) = 1. (7) 
The function tp{x) is the continuous hmit of the vector F defined in section |21 

Remark 1. For equation it can be shown that the former condition is 
automatically satisfied; hence, we can treat it as a problem with no boundary 
conditions. As for equation the degeneracy at the endpoints with such in- 
tegral boundary conditions turn it in a very nonstandard problem in parabolic 
pdes. We discuss some of the issues raised by this problem in section\^ 



4 The Kimura connection 

Using mean-field Gaussian approximations for the frequency independent 
case, Kimura ^Uj has derived a PDE for the evolution of the transient fixation 
probability, which will presumably evolve to a stationary solution that will 
then be the standard fixation probability. This equation is now known as 
the Kimura equation and read as follows: 

dj = x{l-x)dlf + -fx{l-x)dj, /(t,0)=0 and f{t,l) = l. (8) 

The stationary state can be readily found as 



which corresponds to given by (jZj), with a = and /3 = 7. 

Let f = f{x) + fs{x), then / satisfies (jHl) with homogeneous boundary 
conditions. In sectional it will be shown that p can be written as a sum of a 
smooth part q with a distributional part with support at the end points. It 
will be also shown that q satisfies (jH)) without boundary conditions. 

Now let us assume that f{t, x) is sufficiently smooth in x. Then a straight- 
forward computation shows that 

/ [a;(l — x)d1f{t, X + 7a;(l — x)dxf{t, x)~\ q{t, x)(ix = 
Jo 

/ /(t, x) [dl - x)g(t, x)) - 7(9^ - x))] dx. 
Jo 
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Thus, dH} and ^ with no boundary conditions are formally adjoints with 
the appropriate inner product. 

A further relationship between / and q should be pointed out, namely 
that, up to a normalising constant, we have 

fit, x) = x{l — x)q{t, 1 — x). 

The adjointness discussed above also hold when a is nonzero. In this case, 
we have the generalized Kimura equation given by 

dtf = xil - x)dlf + 7x(l - x) (/?+(« - f3)x))dj, 

f{t,0) = and /(t,l) = l. 



5 Mathematical issues 

There are a number of important questions related to equations (0) and 
(0) given the degeneracy at the endpoints and the non-standard boundary 
conditions. 

Note that the last two equations will, generally, have very different qual- 
itative behavior as t — oo. In particular, we prove the following, concerning 
equation (jH)): 

Theorem 1. 1. For a given G L^([0, 1]), there exists a unique solu- 
tion p = p{t, x) to Equation ^ of class C°° x (0, 1)) that satisfies 
p{0, x) = p^{x). 

2. The solution can he written as 

p{t, x) = q{t, x) + a{t)5o + b{t)Si, 

where q G C"^(]R+ x [0,1]) satisfies (0j without boundary conditions, 
and we also have 

(^{t) = / q{s,Q)ds and b{t) = / l)ds. 
Jo Jo 

In particular, we have that p G C°°(M+ x (0, 1)). 
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3. We also have that 

lim q{t,x) = (uniformly), lim a{t) = 7Co[p^] and lim b{t) = 7ri[p°], 

where 7ro[p°] = 1 — 7ri[p°] anc? the fixation probability associated to the 
initial condition is 

lo IyP°ix)dx exp{-y^^ -yP)dy 

Note that this means that the solution solution will 'die out' in the 
interior and only the Dirac masses in the end points will survive. 

4- Write p^ = a^6o + W6i + G L^{[0, 1]) and let Aq be the smallest eigen- 
value of the associated Sturm- Liouville problem (cf. //, we assum 
that g° G L^([0, 1], — x)dx) and if \\.\\2 denotes the corresponding 
norm, then we have that 



|g(t,.)ll2<e-^«*||g°(. 



2- 



Moreover, we always have the following L} bounds: 
(a) 

||g(t,.)lli<e-'»*||g°(.)lli; 



(b) 



(c) 



TX,[p']-e~^-'\\q\.)\\,<a{t)<Mp% 



7ri[/]-e-^«*||g°(.)||i<&W<^ib% 



It is important to note that equation © is not a good long-term approx- 
imation for the discrete process in the case of a Hawk-Dove game, as will see 
that it presents no diffusion to the boundaries. In this case, the final state 
of any non-trivial initial condition will be fully determined by the unique 
non-trivial equilibrium of the game, as the following result shows: 
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Theorem 2. Consider p{t,x) G [L^ (1 H ^)([0, 1]) solution of Equation 
Then 

lim pit, x) = 6^* , 

t—^oo 

where x* = — a) . 
Proof. Consider 

(xia-P)+P)'^ 1 1 a-(3 

Then, x(l — x){x{a — (3) + P)(j)'{x) = —(f>{x), which imphes 

dt / p{t, x)(j){x)dx = — / p{t,x)(j){x)dx , 
Jo Jo 

and we conclude that the final state is supported at x*, the only zero of 4>{x). 
Using the conservation of mass, we prove the theorem. □ 

Remark 2. 1. For the case of non Hawk-Dove game, i.e., a game only 
with trivial stable equilibrium, then we have 

lim p{t, x) = c6o + (1 — c)6i, 

t— >oo 

The constant c is directly related to the fixation probability, in the fol- 
lowing sense. Let vro[po]e be the fixation probability found with a and 13 
replaced by e~^a and e~^P respectively. Then, |^ show that 

c = 7To[p'] + 0{e). 

Thus, if a and (3 in the original problem are interpreted as scaled down 
selection parameters, then equation (0j yields the same asymptotic be- 
haviour. 

2. For initial conditions in -C/"^([0, 1]), an adaptation of the boundary cou- 
pled weak solution developed in may be used to show similar results 
for games with or without a non-trivial equilibrium. 

Equation ® is a good approximation for the discrete case, as can be seen 
in the following: 
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Theorem 3. Let pN^^t{x,t) be the solution of the finite population dynamics 
(of population N, time step At = l/N"^), with initial conditions given by 
p%{x) = X = 0, 1/N,2/N, ■■■ ,1, forp^ e L^([0, 1]). Assume also that 

(A-1,S-1,C-1,L)-1) = l/N{a,h,c,d) + 0{l/N^). Letp{t,x) be the 
solution of equation with initial condition given byp^{x). If we write p^ 
for the i-th component of p]\f^Atix,t) in the n-th iteration, we have, for any 
t* >0, that 

\im p%" =pit,x), xe[0,i], te[o,t*]. 

Equation (0) is however a good approximation of © for intermediate 
times and strong selection. 
In fact, 

Theorem 4. Consider {a', (3') = e{a,l3) and t' = e^^t. Then, in the limit 
e 0, we have that the regular part of the solution of the re-scaled equa- 
tion (0) converges to the solution of equation ^ in L^([0,T] x [0,1], dt ® 
x{l — x)dx), if the initial condition is in H^{[0,l],x{l —x)dx). 



Proof. Dropping and having in mind Theorem^ we re- write Equation © 

as 

dtq, = edl (x(l - x)q,) - (x(l - x){x{a -(3)+ (3)q,) (9) 
and then we have the a priori estimate 



1 

- / x(l — x)qldx 
2 Jo 



e I (9x. (x(l - x)qs)f dx + l- [ (xa + (1 - x)P)d^ ((x(l - x)qef) 



2 Jo 



< — — — [ x{l — x)q'^dx 



We differentiate equation Q with respect to t, proceed as above to find 
the estimate 

1 

x{l-x) (dtqefdx < , 
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for an e- independent function $i. In order to find an e- independent bound 

for x(l — x) (OxQe) dx, first we prove 



Idt [ q^.dx < \ f [x{l - x){x{a - P) + P) - e{l - 2x)] d^qldx <C [ q^^dx 
2 Jo Jo Jo 

and this imphes an a priori bound for q'^dx. Then, note that 

-dt / x{l -x) {d^qefdx < 
^ Jo 



1 r 







-^(a - - x) - (1 - 2x)(xa + (1 - x)/5) 



{d^qef dx 



+ / [(xa + (1 - x)/3) - (a - - 2x)] x(l - x)(9^g^Mx 
Jo 

<Ci [ x(l - x) {d^q.f dx + C2 / g,Mx . 
7o Jo 

We conclude an a priori bound for /q^x(1 — x) (^^..^e)^ dx and then from 
Rellich's theorem, we know that x(l — x) (q'e)^ dx is in a compact set of 
L^([0,T] X [0, 1]). This proves the theorem. □ 

Remark 3. Equation ^ and ^ have a very important difference, even 
in the case where their asymptotic behaviour is the same. For equation ^ 
the Diracs at the endpoints appear at time t = 0"^, while for (0j this is only 
attained at t = 00. Thus, we have the unusual situation that, at the endpoints, 
the parabolic problem is more singular than the hyperbolic associated problem. 

6 The replicator dynamics connection 

The replicator dynamics models the evolution of the fraction of a given type 
of individuals in a infinite population framework. For a pay-ofT matrix given 
by 

;: :) . 

in its simplest form the replicator dynamics reads as following 

X = X{1~ X){X{a- (3) + (3) . (11) 
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Equation ^ can be written as 



dtp + x{l — x){/3 + {a — (3)x)dxP+ 

+ (/3 + 2(a - 2[3)x - 3(a - [3)x^) p = 

Its characteristics are given by 

— =x(l-x)(/3+(a-/3)a;), 
as 

dz 

= -(/? + 2(a - 2/3)a; - 3(a - p)x^) z. 

as 

The projected characteristics in the x x t plane are given by 

dx 

— =x(l-x)(/?+(«-/?)x), 

which is just (fTTj) . 

For smooth solutions, one can then write the solution to (0) — as done in 
[H] — as 



p(t,a;) = a% + 6% + g°($_t(x)) 



P + {a - (3)<l>.t{x) 



/3 + {a- /3)x 



— X) 

(12) 

where is the flow map of (fTH) . 

Notice that, when a — /5 7^ 0, the first order term does not represents a 
pure drift, but also a dampening (enhancing) for a > /3 (a < /5, respectively). 

Thus, equation (0) can be seen as an Eulerian representation of a quantity 
associated to the probability density evolution, but not to the probability 
density itself. If we let q{t, x) = p{t, x) — a^6o — b^6i, we see, from flT^ . that 
the Lagrangian transported by the replicator flow is 

u(t, x) = x(l — x)(/3 + (a — (3)x)q{t, x). 

Thus, (llip can be Langragian representation of m, once the initial 

probability distribution is given. Since, we can recover q from m, and hence 
can recover p, we have that, when there is no diffusion, solutions to ()lip 
together with initial probability distribution are equivalent to (0). 
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An interesting question is to quantify how good is the dynamics given by 
equation (0) — or equation pij) for that matter — as an approximation to the 
dynamics of (jH)) in the case of small diffusion, i.e., strong selection. Besides 
the results already alluded to in section the following results have also 
been shown in f6j: 

1. For games without a non-trivial stable equilibrium, we have that the 
dynamics of p is well approximated by solutions of © over a long time 
modulated by an envelope on a slow timescale. 

2. For games with a non-trivial stable equilibrium, the above holds away 
of such an equilibrium. Near the equilibrium, we have a balance of 
diffusive and selective effects. This prevents the Dirac formation at the 
equilibrium point. 

3. Combining the remarks above, we have, for Hawk-Dove games, that a 
non trivial initial distribution (i.e., that is not peaked at the endpoints) 
tends to peak at the interior equilibrium, and that such a peak takes a 
long time to die out. For an example see figure 11 in [7j 

7 A numerical tour 

For a compairison of the discrete and continuos models, as well as an extensive 
ensemble of simulations for (jH)), the reader is referred to 0. 

Here, we shall focus on compairing the solutions to (jHI) with solutions to 
(0). We present two sets of simulations of (0); with large /3 and large a. We 
then compared the solutions to (0) with rescaled time and coefficients. We 
also plotted the position with the peak — with rescaled height — with the peak 
of the solution to ©• For display conveninence, we have omitted the very 
ends of the interval and plotted Axp instead of p. 

8 Further remarks 

The analogy between the Moran process for finite populations and the repli- 
cator dynamics can be taken further. More precisely, suppose that the indi- 
viduals taking part in the Moran process do not play only pure strategies as 



16 




Figure 2: Solutions to equation ®, labled as Moran, and to (0) labeled as 
Nondiffusive in the frequency independent case, with a = (3 = 20 and initial 
condition p^{x) = x{l — x)/6. 




Figure 3: Evolution of solutions to © together with the peaks given by 
solutions to (0) plotted as points with rescaled height for a convenient display. 
Same paramters and initial condition of figure |21 
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Figure 4: Solutions to equation ®, labled as Moran, and to (0) labeled as 
Nondiffusive in the frequency independent case, with a = —20 and (3 = 20 
and initial condition p^{x) = 20x^(1 — x). 




Figure 5: Solutions to equation ©, labled as Moran, and to © labeled as 
Nondiffusive in the frequency independent case, with a = —20 and (3 = 20 
and initial condition p°(a;) = 20x'^(l — x). 



18 



0.04 
0.035 

0.03 
0.025 

0.02 
0,015 

0.01 
0.005 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Figure 6: Evolution of solutions to © together with the peaks given by 
solutions to (0) plotted as points with rescaled height for a convenient display. 
Same paramters and initial condition of figures ID and 

in the above analysis, but are allowed to play mixed strategies. In particular, 
let us suppose that the game involves two kind of strategists, Eq^ and Eq^, 
where an i^^-strategist means that he/she plays pure strategy I with prob- 
ability 9 and II with probability 1 — 9. Then, the pay-off matrix is given 
by 





Ee^ 


Ee2 


Ee^ 


A 


B 


Ee2 


C 


D 



where 

A := 9lA + 9,{l-9i){B + C) + il-9iyD , 

B := 9i92A + 9i{l-92)B+{l-9i)92C+il-9i){l-92)D , 

C := 9^92A + {l-9i)92B + 9^{l-92)C + il-9^){l-92)D , 

D := 9lA + 92il-92)iB + C) + il-92)''D . 

The associated thermodynamical limit is given by 




dtp = dl - x)p) - 4 {x{l - x){x{9i - 92f{a - 13) + {9i - 02)(^2a + (1 - ^2 
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The final state is given by = 7ro[p°]5o + 7ri[p°]5i, where 7ro[p°] = 1 — 7ri[p°] 
and the fixation probabihty 7ri[]9'^] is given by 

lo F(di,e2){y)(^y 

where 

F^g,^e,){y) ■■= exp (^-y\9, - ^2)'^^ - - ^2)(^2« + (1 - ^2)/3) 

Note that the neutral case (i.e., when the two types of individuals are of the 
same kind) is given by 61 = 62, and in this case the governing equation is 
purely diffusive and fixation probability associated to a given initial state is 
simply given by 

vrf = [ xp\x)dx . (13) 
Jo 

We say that an Eq^ strategist dominates an Eg^ strategist {Eg^ >- Eg J if 
the fixation probability of the first type, for any non-trivial initial condition 
is smaller that the one in neutral case given by equation (fT^. With this 
definition, we can prove the following theorem: 

Theorem 5. Eg^ >- Eg^ if and only if the flow of the replicator dynamics is 
such that 61 — > 62- 

As a simple corollary, we have that ii 6* = [3 / {[3 — a) E (0, 1) (the ESS of 
the game), then Eg* >- Eg, \/9 7^ 9*. This shows that an individual playing 
a mixed strategy with probabilities equal to the one of the game's ESS is 
better equipped to win any context. But, as we saw in the previous sections, 
for populations of pure strategists we can not expect an stable mixture (even 
in fractions equivalent to the game's ESS) to evolve. 
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